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1 The sum rules for Fp^. Theoretical Status. 

I start with the consideration of the sum rules for the first moments of the spin structure functions 

1 

rp,n(Q') = / dxg,p,n{x,Q^) (1) 



I will discuss the uncertainties in the theoretical predictions for Tp^n and compare the theoret- 
ical expectations with experimental data. The aim of this consideration is to obtain from the 
experiment the restrictions on the uncertainties in the theoretical description of the problem. 

Consider first the Bjorken sum rule [1], which now is a corner stone not only of the prob- 
lem of nucleon spin structiu^c functions, but of the whole theory of deep inelastic lepton-hadron 
scattering in QCD. (It is interesting to note that in' his original paper 1^1 in 1966 Bjorken had 
classified this sum rule as a " worthless equation" ) . 

The Bjorken sume rule reads: 

TAQ') - Km = 1^41 - ^ - 3.6 (^)' - 20 (^)'l + (2) 

where is the axial /9-decay coupling constant and the last term in (2) represents the twist-4 
correction. The perturbation QCD corrections are known up to the third order '^'^^ (there is also 
an estimate of the fourth order term [^], which is not included in (2), since the uncertainties in 
the included terms are larger than its contribution). The coefficients in (2) correspond to the 
number of flavours Nf = 3. (In the domain of existing experiments only three flavours of quarks 
are effective). In what follows 1 will so often transfer the data to the standard reference point 
Qq = lO.bGeV'^ - the mean value of Q^, at which EMC and SMC experiments were done. 

Let us discuss the perturbative corrections in (2). Today there is a serious discrepancy in the 
values of found from different experiments. The average value of a^, obtained at LEP is [^] 

a^(m|) = 0.124 ±0.007 (3) 

In two-loop approximation this value corresponds to the QCD parameter A3 for three flavours 

A3 = 430 ± lOOMey (4) 

On the other hand,Tte data on the T — > hadrons decay give '^1 (the first ag correction '^1 is 
accounted) 

Q;^(m^) = 0.178 ±0.010 (5) 

from which it follows 

A3 = 170 ± SOMeV (6) 

The small error in (6) is caused by the fact that the partial width r(T 3g), from which (5) 
was determined, is proportional to and the ag correction to it is small. A strong contradiction 
of (4) and (6) is evident. 

The overall fit 1^1 of the data of deep inelastic lepton-nucleon scattering gives in the NLO 
approximation A3 = 250MeV (the error is not given). New data in the domain m| indicate 
lower «,(m|) (SLD I^l : 0.118 ± 0.013,0.112 ± 0.004; OPAL I^l; 0.113 ± 0.012), but the data of 
AMY on e+e- annihUation at ^/s = 57.3GeV results in as{m%) = 0.120 ± 0.005. Finally 
from r- decay it was obtained t^^^: 

a,(m^) = 0.33 ±0.03 (7) 
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corresponding to 

A3 = 380 ± GOMeV 

But the determination of ag from r-decay can be criticized '^^1 on the grounds that at such a low 
scale exponential terms in may persist in the domain of positive — rn^ besides the standard 
power-like terms in accounted in the calculation. 

In such a confusing situation 1 will consider two options - of small and large perturbative 
corrections. In the first case I will take A3 = 2QQMeV . Then 

as{Ql) = 0.180 ± 0.010 Tp{Ql) - Tn{Ql) = 0.194 (8) 

In the second one A3 = AQQMeV and 

a,{Ql) = 0.242 ± 0.025 V^) - -^n{Ql) = 0.187 (9) 

The twist-4 contribution was disregarded in (8) and (9). There is also a discrepancy in its 
value. The value of 6p_„ was determined in the QCD sum rule approach by Balitsky, Braun and 
Koleshichenko (BBK) l^^]. 

Vn = -OMbGeV'^ (10) 

On the other hand, the model [i"^^!*^! based on connection ^^''1 of the Bjorken sum rule at large 
with the Gerasimov, Drcll, Hcarn sum rule at = '^^1 gives 

= -0.15GeV^ (11) 

In what follows I will consider (10) and (11) as two options which correspond to small (S) 
and large (L) twist-4 corrections. I will discuss both approaches of determination of higher 
twist corrections in more details below. From my point of view, no one of these approaches is 
completely reliable and I will use them in comparison with experiments only as reference points. 



I turn now to the sum rules for Tp and F^: 



^_ 3.6(^)2 -20(^)3 

TT TT 71 . 



{±ga + ^as)+ (12) 



f -0.55(^) 

671 TT 



(The correction to the singlet part was calculated in '^^1). 

According to the current algebra qa, as and E are determined by the proton (or neutron) 
matrix elements of flavour octet and singlet axial currents 

- 2msu,a8 =< P, s | j^J \p,s> - 2ms^E =< p, s | ji^J \p,s>, (13) 



-2ms A =<P,s\ f^'J \p,s> (13') 



where 



In the parton model qa, dg and E are equal to 

qa^ ^u- Ad as^ Au + Ad- 2AS E = + Ad + As, (14) 
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where 



A,^/b.W-,-Wl ,^M. (15) 



and q± are the quark distributions with the spin parallel (antiparallel) to the proton spin, which 
is supposed to be longitudinal (along the beam). Ag in (12) has the similar meaning, but for 
gluons, as Ag for quarks 
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Ag{Q') ^ jdx [g+{x, Q') - g^x, Q')] (16) 





Unlike gA, Os a-nd E which in the approximation used above are independent and have zero 
anomalous dimensions, Ag anomalous dimension is equal to -1. This means that 

^9{Q')q^^oo ^ clnQ^ (17) 
The conservation of the projection of the angular momentum can be written as 

^E + A^(g2) + L,(g2) = l (18) 

where has the meaning of the orbital momentum of quarks and gluons. As follows from 
(17), (18) at high Q"^ Lz{Q'^) must compensate Ag{Q'^), Lz{Q'^) ~ —clnQ'^. This means that the 
quark model, where all quarks are in jS'-state, failed with increasing. 

The gluonic contribution to Tp^n term, proportional to A (7 in (12), was calculated in [2-20-23] _ 
There was a wide discussion in the past years if gluonic contribution AgTp^n to Tp^n is uniquely 
defined theoretically or is not [23-36] ^^^q problem is that gluonic contribution to the structure 
functions, described by imaginary part of the forward 7^irt~gluon sccattering amplitude (Fig.l) 
is infrared dependent. Since in the infrared domain the gluonic and sea quark distributions 



Fig.l 

The photon-gluon scattering 
diagrams, the wavy, dashed 
and solid lines correspond to 
virtual photons, gluons and quarks. 



are mixed and their separation depends on the infrared regularization scheme, a suspicion 
arises that A^Fp „ can have any value. This suspicion is supported by the fact that in the lowest 
order in ag the terms proportional to InQ^ are absent in AgVp ^ and this contribution looks 
like next to leading terms in nonpolarized structure functions where such an uncertainty is well 
known. 

In the framework of the operator product expansion (OPE) there is only one operator - the 
singlet axial current jj^^ - corresponding to the first moment of gip + gin, i-e., Fp+„. This fact 
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can be also used as an argument in the favour that (in the ag and higher orders) the separation 
of the terms proportional to E and A in (12) is arbitrary. 

In order to discuss the problem consider the gluonic contribution gip{x, Q^)gi to the proton 
structure function gip{x,Q^), described by the evolution equation 




9ip{X:Q%^Nf^— / ^A[- [g4y,Q')-g_{y,Q% (19) 



where the asymmetry A{xi) is determined by the diagrams of Fig.l. The calculation of the 
asymmetry A{xi),xi = —q^/2pq results in appearance of integrals 



/ fhl , ri=l,2, 

J {k^ — xi{l — xi)p'^ + rrig)"'^ ' ' 



which are infrared dependent. To overcome this problem it is necessary to introduce the infrared 
cut-off (or infrared rcgularization), to separate the domain of large k'j_, where pcrturbativc QCD 
is reliable, from the domain of small k"^. The contribution of the latter must be addressed to 
noncalculable in perturbative QCD parton distribution. Such a procedure is legitimate because 
of the factorization theorem which states that the virtual photoabsorption cross section on the 
hadronic target h, cr1{x, Q^) can be written down in the convolution form 

al{x,Q')=J2<^n^,Q',M')®Mx,M'), (20) 

i 

where a] is the photoproduction cross section on the i*'* parton {i = q,q, g), fi/h are the parton 
distributions in a hadron h, <S> stands for convolution. Both aj and fi/h depend on the infrared 
cut-off M^, but the physical cross section a^/i is cut-off independent. The variation of 
corresponds to redistribution among partons: the trade of gluon for sea quarks. 
As follows from (19) 

< > 



A,rp^Nf^^AgA{M^), (21) 
where 

A(M^) = [' dxiA(xi, M^) (22) 
The convenient way is to introduce cut off in k^^ 

kl>M\x,,p') (23) 

Generally, may depend on x\ and p^. For example, the cut-off in quark virtuality in Fig.l 
-k"^ > Ml corresponds to the form (23) with = (1 - xi){Ml +p'^Xi) if Xi < -Ml/p^. 

In the calculation of the diagrams of Fig.l it is reasonable to neglect the light quark masses 
in comparison with gluon virtuality p^, since we expect that | | is of order of characteristic 
hadronic masses, | |~ IGeV"^ l^^l. Then, introducing the infrared cut-off (23) we have '^^1 



A{M^) = i^l - J dxi{l - 2xi)[lnr - r]^ 



(24) 



where 

xi(l-a;i)p2-M2(xi,p2)- ^ ' 
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If = Const - a rectangular cut-off in k\- , the integral in (24) vanishes (the integrand is 
antisymmetric under substitution xi — > 1 — xi. Then A = —as/^n and we obtain the gluonic 
contribution to ^p,n (12). However, other forms of M^(xi,p^) result in different values of ^(M^), 
what supports the claim l^^l that A[M'^) is cut-off dependent. Even more, if we put = 0, 7^ 
- the standard regularization scheme in the calculation of nonpolarized deep inelastic scattering - 
we will find A — O'^^l. This result is, however, nonphysical because the compensation of —as/27r 
term in A arises from soft non-perturbative domain of /c^ ~ m^, which must be attributed to 
sea quark distribution. 

Although generally A, as well as Ag and AgTp are infrared cut-off dependent, a special class 
of preferable cut-off's can be chosen. In OPE the mean asymmetry A is proportional to the 
one-gluon matrix element of axial current 

r^Aa(0,p,p) =< g,ex I i/.5(0) \ g,ea> (26) 

at zero momentum transfer. It can be shown '^^1 that this quantity is proportional to the 
divergence /^r^Acr(^Pi,P2), / = Pi — P2 in the limit —>■ 0, i.e. to the one gluon matrix element 
of axial anomaly. 

The proof is the following. The general expression of the matrix element of axial current over 
gluonic states with nonequal momenta has the form '^^1 



^ (27) 



^nX<jp\Pl +P2)p + -^^tiapvP\pP2u — -^^ix\pvPlpP2v 



where it was assumed that Pi = P2 = p^ < 0. When deriving (27) we used only Lorenz 
invariance, Bose symmetry of gluons and the conservation of vector currents at gluonic vertices. 
The formfactors Fi_2(^^,P^) have no kinematical singularities if 7^ and, particularly, Fi(/^,p^) 
has no pole at l"^ — 0. Prom (27) we have: 



^Xcpr Pip P2t 



r^xa{0,p,p) = F2(0,p')e^Aap Pp (28) 

and in the limit — the mentioned above statement follows. 

Let us recall the derivation of the anomaly in perturbation theory t^^l . The anomaly corre- 
sponds to the diagrams of Pig. 2. 



Pig.2 

The diagrams for the gluonic matrix 
element of axial current. 



These diagrams have a superficial linear divergences. Cancellation of the divergences in the 
diagrams Pig's 2a and 2b requires shifting of integration variable in 2b: 

k-pi-p2 
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This shift corresponds to substitution xi ^ 1 — xi in the integrated over xi diagrams of Fig.l, i.e. 
in eq.(24). We came to the conclusion that in order to retain connection of r^;^o.(0,p,p) with the 
anomaly which follows from (28) and to preserve the standard form of the anomaly, the infrared 
cut-off M'^{xi,p'^) must be symmetric under Xi ^ 1 — Xi. In this case A = —as/2'7r and we 
come to the gluonic contribution to Tp^n given by eq.l2. With such a cut-off the anomaly comes 
entirely from the hard momentum region and can be considered as a local probe of gluon helicity. 
(It must be mentioned, however, that any other forms of infrared cut-off are also possible, but 
they are less attractive). Some care, however, is necessary, when the calculations with this 
cut-off are compared with HO calculation in nonpolarized scattering, where as a rule, another 
regularization procedure —p^ = 0, 7^ is used. 

The gluonic contribution to the spin dependent proton structure function gip{x, Q^)gi can be 
found experimentally by measuring inclusive two jets production in polarized DIS ^^1. However, 
only large —k"^ component of Fig.l diagrams can be determined in this way, the small — fc^^ 
component is nonmeasurable, since in this case it is impossible to separate one jet from two 
jets events '^^l This fact is in complete accord with the formulated above statement about the 
arbitrariness of infrared cut-off. 

Now about the numerical values of the constants gA, os and E entering eq.(12). is known 
with a very good accuracy, gA = 1.257 ± 0.003 '^l.It must be mentioned, however, that eq.(13') 
follows from the assumption of exact isospin symmetry. In fact, one may expect its violation of 
the order of 1%. Therefore, the uncertainty in this number is not given by experimental error, 
but may be ~ 1%. Under assumption of SU(3) flavour symmetry in baryon decays ag is equal 
to 

ag = 3F - D = 0.59 ± 0.02 (29) 

where F and D are axial coupling constants of baryon /?-dccays in SU(3) symmetry and the 
numerical value in the r.h.s. of (19) follows from the best fit to the data t'^^'. The combination 
3F — D can be found also from any pair of baryon /3-decays. The comparison of the values of 
3F — D, obtained in this way shows, that the spread is rather narrow \ Sag \< 0.05 and at 
least as > 0.50. This may be an argument in the favour that SU(3) violation is not large here. 
(For the recent suggestion of such violation see, however, t'^'^l). It should be mentioned that at 
fixed Tp n the uncertainty in ag only slightly influences the most interesting quantities E and As. 
As follows from (12) and (14) 

(5E = -S{As) = ^Sas. (30) 

as was also determined by the QCD sum rule method '^^l In this approach no SU(3) flavour 
symmetry was assumed and the result 

= 0.5 ±0.2 (31) 

is in agreement with (29), although the error is large. 

What can be said theoretically about S? In their famous paper l^^l Ellis and Jaffe assumed 
that the strange sea in the nucleon is nonpolarized. As = 0. Then 

E^as^ 0.60 (32) 

(The sum rule (12) in the framework of this assumption is called Ellis- Jaffe sum rule). This 
number is in a contradiction with the experimental data pointing to smaller values of E. On 
the other hand, Brodsky, EUis and Karliner I^^l had demonstrated that in the Skyrme model at 
large number of colours Nc, E ~ 1/A^c and is small. From my point of view this argument is not 
very convincing: the Skyrme model may be a good model for description of nucleon periphery, 
but not for the internal part of the nucleon determining the value of E (see also ^^^^). 
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2 Calculations of the Matrix Elements over the 
Polarized Nucleon by QCD Sum Rule Approach 



The QCD sum rule method was used to calculate the nucleon coupling constants gA'\'^s^^ 
the coefficients hp^n which determine the twist-4 contributions I^^l. The basic features of the 
approach are the same as originally suggested for calculation of nucleon magnetic moments 

[46,47] _ 

We wish to find the diagonal matrix element 

<p,s\ i(0) \p,s> (33) 

To this end add to the QCD Lagrangian the term 

AL = i(x)5, (34) 

where 5 is a constant external source. In the case of qa determination for j and S we substitute 
j^5^ and A^^^ - the constant external axial field etc. Consider the polarization operator of the 
currents r] with the nucleon quantum numbers 

n(p) = i J d'^x e'P'' < I T {r]{x), f]{0)} \ > (35) 

rjp = {u''{x)Cjy{x)) 7^75(^'(2:)e"'^ (36) 
Vn = Vp{u ^ d) 

where a, 6, c = 1, 2, 3 are colour indeces. Separate in n(p) the term proportional to external 
source S: 

r{p,p,0)S^f j d^xd^ze'^'^S <\ T{r]{x),j{z),fi{^)} \ > (37) 

Suppose that < 0, | |S> where Rc is the confinement radius, and perform OPE in 
Unhke the standard OPE for polarization operator, the presence of external source results 
in appearance of a new type of vacuum expectation values (v.e.v.), induced by external source. 
If the source 5" is an axial field A^, then, e.g. 

< I M7^75M I >A= cA^ (38) 

It can be easily shown that for massless quarks iriq = {q = u,d, s), c — in the case of 
isospin vector or octet A^ axial field. In the first case we can write 

1 f 

< I M7^75M I >A= fim - J d^xe'*^ x 

X < I T {jS(x), jg)(0)} I >= limng(g)4^) (39) 
Since the axial current jj^^ conserves 

nS(g) = -(W-9/.9.)n(g^) (40) 
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The non-zero result for n^j,(g) at g ^ comes from the pole ~ in n(g^) which corresponds 
to the massless pion intermediate state 

ng(^)=(-^ + 5..)/.^ (41) 

In (41) the pion mass in the propagator was accounted as the most important effect of nonva- 
nishing quark masses at small Substituting (41) into (39) and going to the hmit — > we 
obtain the desired result. 

The same proof may be repeated for J^^^ in the limit of SJJ{?>) flavour symmetry. The 
vacuum expectation values which enter the OPE of the r.h.s. of (36) can be classified according 
to their dimensions and the number of loops (the dimension of = ^Al^^ is 1, the chirality 
conserving structure is considered) 

Dimension V.e.v. Number of loops 

1 \.A^ 2 

3 < I ^7^75? I >A= 1 

(42) 

5 < I a,G\^ I > 2 

5 < I g7.G'IJ4A«g I >A= f^A^ 1 

7 < I I >2 

Here G^^ = {l/2)e^^xaGxa, rn\ ^ 0.2 GeV^ t^^'. The important role of induced by external 
field v.e.v. is evident - besides the unit operator they correspond to lowest dimensions and to 
minimal number of loops. The OPE for the r.h.s. of (37) can be constructed exploiting the v.e.v. 
given in (42). 

In order to find the matrix element (33), represent the l.h.s. of (37) in terms of contributions 
of physical states using dispersion relations. Generally, when the momentum of external field 
is non-zero, P is a function of three variables, P = V^p^^p^] q^), and may be represented by the 
double dispersion relation: 

oo oo / 2\ 

nplphq')- [ ds, j ds, . P^'^-J'^^ ,, +Pipl)f{pW)+ (43) 

^ ^ l-^l - Pl)\-'^2 - P2) 

+p{pI) fiPi e) 

, where P{p^) is a polynomial and /(p^, ?^)is given by the ordinary dispersion relation 

^ — — + subtr. terms (44) 

s — 
^ 

We are interested in the limit q ^ O-pf = pi = and at the first sight it seems that one 
vartiable dispersion relation for T{p'^,p'^, 0) can be written in this case. Indeed, 

7 ,1,, J is, , = / ^<£ii^d.,ds, ( ^ L_) (45) 

J J [Sl - P^)[S2 - P^) J S1-S2 \S2-P^ SI-PV 
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In the first (second) term in the r.h.s. of (45) integration over Si(s2) can be performed and 
the result has the form of one-variable dispersion relation. Such transformation is, however, 
misleading, because, in general, the integrals 



dsi 



P{S1,S2,0) 




dS2 



p(si,S2,0) 



Si - S2 



Si - S2 



are ultraviolate divergent. This ultraviolate divergence cannot be cured by subtractions in one- 
variable dispersion relation: only subtractions in the double dispersion representation (43) can 
be used. It is evident that the procedures which kill the subtraction terms and lead to fast 
convergence of dispersion integrals in standard one-variable dispersion representations like Borel 
transformation in do not help here. 

Let us duscuss the determination of the coupling constant as in more details. Consider in 
(37) the coefficient function at the structure {A^^^p^)p^5. The advantage of this structure is that 
it has the largest power of momenta in the numerator. Therefore, the integrals for the coefficient 
function converge better and the contribution of excited states in the dispersion relation (43), 
which is the background in our calculation, is smaller. It can be easily seen that the function p 
in the r.h.s. of (43) which correspond to the bare loop diagram Fig. 3 has in this case the form 



where a is a calculable constant 



Fig. 3. The bare loop diagram, corresponding to determination of the couphng constants or ag. The sohd 
lines correspond to quark propagators, crosses mean the action of currents t], f), the bubble corresponds to quark 

interaction with external field. 

The substitution of (46) into (43) gives for the first term in the OPE in the r.h.s. of (37) at 



In this simple example the dispersion representation is reduced to one- variable dispersion rela- 
tion, but with the square of (si —p^) in the denominator. Of course, by integrating by parts (47) 
may be transformed to the standard dispersion representation. However, the boundary term 
arising at such transformation must be accounted; it does not vanish even after application of 
the Borel transformation. This means that even in this simplest case the representation (43) is 
not equivalent to one- variable dispersion relation. 



p(si,S2) = asiS25{si - S2), 



(46) 



222 
Pi = P2 = P 




,2 



(47) 
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Fig.4. The integration domains the in si, S2 plane 



Let US represent r(p^,p^; 0) in terms of contributions of hadronic states using (43) and separating 
the contribution of the lowest hadronic state in the channels with momentum p I^^l . As is seen 
from Fig.4, it is convenient to divide the whole integration region in si,S2 into three domains: 
I) < si < 1^^ < S2 < W^; II) < si < W\ < S2 < 00; < si < 00, < S2 < 

W^; I1I)W^ < si < 00, <S2<oo. Adopt the standard in QCD sum rule model of hadronic 
spectrum: the lowest hadronic state nucleon plus continuum, starting from some threshold W'^. 
Then in the domain I only the lowest hadronic state N contributes and 



where m is the nucleon mass A is the transition constant of the nucleon in the current 77: 
< N \ fj \ >= Xv where v is the nucleon spinor. In the domain III the higher order terms in 
OPE may be neglected and the contribution of hadronic states is with a good accuracy equal to 
the contribution of the bare quark loop (like Fig.3) with perturbative corrections. The further 
application of the Borel transformation in essentially suppresses this contribution. 

The consideration of the domain II contribution is the most troublesome and requires an 
additional hypothesis. Assume, using the duality arguments, that in this domain also, the 
contribution of hadronic states is approximately equal to the contribution of the bare quark 
loop. The accuracy of this approximation may be improved by subtraction from each strip of 
the domain II of the lowest hadronic state contributions proportional to 5(si — m^) or 5{s2 — m^)- 
The terms of the latter type also persist in the functions f{Pi), fiPo) in (43). They correspond to 
the process when the current fj produces the nucleon N from the vacuum and under the action 
of the external current j the transition to excited state N — >■ N* occurs or vice versa (Fig. 5). 



Fig. 5. The schematical representation oi N ^ N* {N* —* N) transitions in the external field. 
At Pi = pi = p^ these contributions have the form 



^(-51, S2) = asX'^6{si - m^)6{s2 - m^) 



(48) 




(49) 
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with some unknown function b{s). The term (49) will be accounted separately in the l.h.s. of 
(43). I stress that the term (49) must be added to the l.h.s. of (43) independently of the form 
of the bare loop contribution p(si,S2). Even if p(si,S2) — 0, when the OPE for the vertex 
function r(p^,p^, 0) with zero momentum transfer starts from condensate terms - the term (43) 
may persist. (49) may be written as 

/■°° /I 1 \ 1 

/ dsb{s) + (50) 

The functions /(p^) in (43) can be represented by dispersion relation as 

m - r -^.ds (51) 

JO s — p'^ 

The integration domain (51) may be also divided into parts < s < W"^ and W"^ < s < oo. 
According to our model the contribution of the first part is approximated by A^-state contribu- 
tion, the second one by continuum. These two parts look like the contributions of the first and 
the second terms in the bracket in (50). 

Now we can formulate the recipe how the sum rule can be written. At the phenomenological 
side - the l.h.s. of the sum rule - there is a contribution of the lowest hadronic state A^ and 
the unknown term (50), corresponding to nondiagonal transition A^ N* in the presence of 
external field; 

{p^ — w?Y 

The contribution of continuum corresponding to the bare loop (or also to the higher order terms 
in OPE, if their discontinuity does not vanish at s — oo) is transferred to the r.h.s. o the sum 
rule. Here it is cancelled by the bare loop contribution from the same domain of integration. 
As a result, in the double dispersion representation of the bare loop the domain of integration 
over si, S2 is restricted to < si, S2 < W^. Finally, apply the Borel transformation in p^ to both 
sides of the sum rule. In the r.h.s. - QCD side - the contribution of the bare loop has the form 



r,,6w^M-^+^ (52) 

Jw^ s — \ p — m'^ s — 



JO Jo Si — S2 

= 2P r'ds^r ds/-!^e-"'-' (53) 

Jo Jo Si — S2 

where P means the principal value and the symmetry of p{si,S2) was used. The l.h.s. of the 
sum rule is equal to 



A' 



«8T7?e-"^'^'^'-^e-"^'/'^'+e-"^'/^' T dsb{s)^^exp[-{s - m^)/M'] (54) 
' Jw^ s — 



where 

r ds-^ (55) 
JvV^ s — 

In (54) A is an unknown constant which can be determined from the same sum rule exploiting 
the fact that the dependence of the second term in (54) differs from the first. Of course, 
the calculation is reliable only if the contribution of the second term in (54) is smaller, say, less 
than 30% comparing with the first. The last term in (54) is of the same origin as the second 



11 



one - it comes from inelastic transitions N ^ N* (Fig. 5). but it is exponentially suppressed in 
comparing with the second and as a rule may be neglected. 

Omitting the details of the calculation, I present the result - the sum rule for determination 
of the coupling constant Og [^^^ 



M2 j 

) L^l' + f-'f.M^E, L-/^ + f (2.) V>?^o (^) L-/^ (56) 

a = -(27r)^ < I I >= 0.55GeV^, b = (27r)^ < | —G^ \ >= 0.5GeV^ 



0^/4 8 
+ ^ -- + 



M2 V 3 9 

Here 



Eo{x) = 1 - e-^ £;i(x) = 1 - (1 + x)e-^ E2{x) ^ 1 - + x + — j (57) 

'"(f) /'"(f) 

For the nucleon coupling constant = 327r'^A^ and continuum threshold we took the values 
found in the mass sum rules [sO'^e]^ _ 2.1GeV^,W'^ — 2.3GeV'^. The normahzation point // 
was chosen as /i — O.bGeV. The sum rule (56) may be compared with the sum rule determining 
the nucleon mass I^^'^^l: 



X%^,e--''"' . M^E. (^) L-'" . 1 < I I > (^) ^ i^^V, 

Adding (56), (59) and applying to the sum the differential operator 

(1 _ M'^)M'e-^'/^' 
which kills the constant A, we get the sum rule for ag: 



(59) 



a8 = -l + -i(l-M^-^ 



^U^'flmlM'E, L-«/9 + a'hA (60) 



The dependence of the r.h.s. of (60) is plotted in Fig. 6. 



Fig.6. 

The Borel parameter dependence of the r.h.s. of (60). 
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It is seen that the dependence in the interval 0.8 < < l.SGeV'^, where the nucleon 
mass sum rule holds, is remarkable. Also the constant A found from (56) is not small enough, 
A ~ — 0.3. For these reasons the error in the determination of ag is large and the final result is 

as = O.btrd (61) 

I would like to mention that the uncertainty in the value of ag arising from proportional to A term 
in the l.h.s. of (56) manifests itself as well in the lattice calculations of these or similar quantities. 
In the lattice calculations the polarization operator (35) is measured at large euclidean times t 
and the result is proportionl to 

~Xl e-""^ [as{l + mt) + 2Am^] (62) 

(see the second reference in t^^l where a similar formula was obtained in the case of nucleon 
magnetic moments). The last term is usually neglected in lattice calculations. As is seen from 
(62) and from the values of ag and A presented above, in order to find ag with 10% accuracy it 
is necessary to go to t > 10m~^ ~ 2/m, where the whole effect is very small. 

I dwell now on the attempt to determine the proton singlet axial coupling constant Qq — 
by QCD sum rules '^^1. The difference in comparison with qa and ag determination comes from 
nonconservation of singlet axial current caused by the anomaly 

d,3^S = ^ + 2^m, s^,s (63) 

In the r.h.s. the strange quark mass term is accounted, whose contribution will be essential. The 
flavour singlet pseudoscalar meson [rj') is not a Goldstone and the used above method (eqs.(38)- 
(41)) of determination of induced by the axial field v.c.v.'s fails. So, a special investigation is 
needed in order to find the value of quark condensate induced by axial field. Like (39) we can 
put 

< I «7^75« I >A= limg^o n(.°j(g)4°) (64) 

where U^^^ is defined by eq.39 with substituting jj^^ (2/3) j'^^ . But now 11^^^ contains both - 
transverse and longitudinal parts 

nW(g) = -UL{qy,. - nt{q'){S,,q' - q,q, (65) 

Uiiq^) and Ilt{q^) have no poles at = and the interesting for us quantity is 

/o' = -n^(0) = [Ul{0) - UL,pertm (66) 

The perturbative part was subtracted in (66) since it will be accounted in explicit way by 
perturbative calcuUation (practically, it is small). The constant /q plays here the same role as 
in (38). 

In order to separate ni(g^) multiply Il^i^ by qfj,qu- We have 

g^g,nW(g) = -q'ULiq') =i^j d'x e^'^^ < \ T {g';,,{x)g;,{x) , 

Glfsmif,{0) + 2tms 5755(0;)} I > (67) 



3a;., 



47r 
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Represent U^^^ (q) in terms of contributions of physical states using the dispersion relations. Then 
it can be easily seen that UlIq^) is contributed by pseudoscalars, the axial mesons contribute to 
Ilt{q^)- Using the notations 

< I q^u^^q \v' >= Wl' (Iv < I \r]' >= ^fn'^l^ (68) 
{q — u, d) and representing ni(g^) in terms of r]' contribution and continuum, we write 

' so 

where sq is the continuum threshold. In (69) (5{s) is determined by perturbative calculation of 
the bare loop (Fig.7) corresponding to the r.h.s of (67). 



Fig. 7. The bare loop contribution to the r.h.s. of (67), 
the dashed Hnes are gluons. 



Therefore, IlL{q^) is given by 

9 so 

m^, — q^ I s' — q" 



n.(5^) ^ n.(5^) - n. ^^(f ) = - / (70) 



To find the quantity consider the integral over the contour C in Fig. 8 in the complex 

plane q^: 



Fig. 8. The integration contour in the complex plane in eq.(71). 



_L|n,(,.)^.n,(0)._,;/,.4/^wf (71) 



C 

The calculation of the diagram Fig. 7 gives 



3 
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To find gl^j/fn' let us write OPE for the r.h.s. of eq.(67) and use the QCD sum rule method. In 
OPE take into account, besides the bare loop, the contribution of the gluon condensate and of 
the quark-gluon condensate (Fig. 9) 

- ^ < I saa/sly G^^ s I >= < I ss I >, (73) 

ttIq = Q.SGeV^, stemming from the term, proportional to in (67). 



Fig. 9. Gluon (a) and quark-gluon (b) condensate contributions to the r.h.s. of eq.(67). 



Let US also account the condensate ~< | g^G^ \ >, taking its estimate from Ref.52. 
The sum rule for g'^//,,' is given by: 



e \ 4 

Here e ^ 0.2GeV^ I^^l 



a, = -{2nf < I ss I >Ri 0.8a ^ 0.44^6^^ (75) 
From (74) it numerically follows 

y^, /,, = (4±l)-10-^Gey^ (76) 

and 

= -UUO) = g;. U' - ^ So « 3.5.10-=^Gey2 (77) 

in comparison with the similar constant for the octet current = (0.133 GeVY = 18.10"^ GeV^. 

Nevertheless that all looks okay in this calculation, the result is wrong. In order to demon- 
strate this, determine the s-quark coupling constant g^, f^i. The only diference with the previous 
calculation will be that in the r.h.s. of eq.(67) the equal-time commutator, proportional to russs 
will appear and, as a consequence, owing to the s-quark condensate an additional term appears 
in the r.h.s. of the sum rule (74) 

1 gm^,/A^^4 ^ Q I I Q g_^Q-3 Q^y2 ^^g^. 
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Thus, our calculation results in a wrong ratio of 77' interaction constants with u (or d) and 
s-quarks 

4 ~ 2.5 (79) 
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whereas it is known that 77' is a flavour singlet and this ratio must be close to 1. 

We arrive at the conclusion that the sum rule for the longitudinal polarization operator with 
the singlet axial current does not work at the values of the Borel parameter of order of 
the Vj mass square. The only way to avoid this discrepancy is to assume that in the OPE in 
this case there are some important higher order terms. This conclusion is not surprising:it has 
been known long ago I^^'^^l that with the standard OPE it is impossible to describe the Okubo- 
Zweig-Iizuka rule violation in the pseudoscalar and longitudinal axial channels. It is necessary 
to take into account higher order terms and instantons in the direct channel were proposed as 
possible candidates for such terms. Now this idea is strongly supported by the calculations in 
the instanton liquid model ^^^1, as well as by the lattice calculations 

Since the attempt to calculate the induced by singlet axial vacuum condensate fails, it is 
impossible to flnd by QCD sum rules the coupling constant a^ — Yu- the part of the proton spin 
projection carried by quarks. In fact, the situation is even wosc. It can be shown '^^1 that if 
the value of S is taken from experiment, then the sum rule or its determination (like (60) with 
the induced by the field v.e.v. considered as free parameters) and the sum rule for transverse 
singlet axial current polarization operator nt(g^) are in contradiction . This indicates the OPE 
breaking down at virtualities ~ IGeV^ in the vertex function for proton interaction with singlet 
axial current or/and in the transverse singlet axial current polarization operator. Thus, the 
sitTiation in the singlet axial channel resembles the one in the pseudoscalar channel and one may 
expect a noticeable violation of the OZI rule here too. A similar trouble, perhaps, faces attempts 
to determine E using the so called "[/(I) Goldberger-Treiman relation" (for a review see l^^l). 
So, at this stage, the only way to find E is from experiment, exploiting eq.l2. 



3 Twist-4 Corrections to Fp^ from QCD Sum Rules 

The general theory of twist 4 corrections to deep inelastic scattering on polarized nucleons has 
been developed by Shuryak and Vainstein t^^l. The have found (see Ref.l3, Errata for corrections 
of errors) 



(r -\-r \ - -- - 

p ^ njtwisti g 



^9 02 J ^ 9l,p±n{x) 



m 



(80) 



Q 

u 

Here indeces S,NS correspond to + and - signs in the l.h.s. of (80), — 5/18, C^^ — 1/6. 
The reduced matrix elements {{U)),{{V)) are related to the matrix elements of the operators 



in the following way 



= ug Gl^i.-X^u 



{N\U,\N)^ s,{{U)) 



(81) 
(82) 

(83) 
(84) 
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where Sfj_ is the unit nucleon spin vector, A^^iy and Si,,a stand for (anti)symmetrization over the 
given subscripts. The indeces S, NS mean 

18 - 
S ^ UU + dd + —ss, NS uu — dd (85) 
5 

{{U)),{{V)) were calculated by Balitsky, Braun and Kolesnichenko (BBK) l^^J using the QCD 
sum rule approach. The result for was given in (10). For it was obtained 

bp+^ = -0.022 GeV^ (86) 

This result, however, cannot be considered as rehable for the following reasons: 

1. BBK use the same hypothesis as EUis and Jaffe did, i.e., assume that s-quarks do not 
contribute to the spin structue functions and instead of singlet (in flavour) operator consider the 
octet one. 

2. When determining the induced by external field vacuum condensates, which are very 
impotant in such calculations they saturate the corresponding sum rule by r] meson contribution, 
what is wrong. (Even the saturation by ry'- meson would not be correct since rj' is not a Goldstone) . 

3. One may expect that in the same way as in the calculation of S by the QCD sum rule, in 
this problem the OPE series diverge at the scale ~ 1 GeV where the BBK calculation proceeds. 

Even in the case of the Bjorken sum rule, where the mentioned above problems are absent, 
the value 6p_„ obtained by BBK is questionable. Let us consider this calculation in more details. 
The bare loop diagram for this case is shown in Fig. 10. 



Fig. 10. The bare loop diagram for twist 4 correction to the Bjorken sum rule 
for deep inelastic electron-nucleon scattering. The dashed line corresponds 
to discontinuity over pi at pi ^ P2- 



In [^^1 this diagram was calculated by introducing an ultraviolet cut off n"^. It was found 
that for the chosen Lorentz structure the singular in p'^ term is proportional to ln^{fi'^ — p^). 
Such cut off dependence reflects the fact that the spectral function p(si,S2) in eq.(43) is not 
proportional to 5-function. The logarithm square dependence of T^p^) on the cut-off cannot 
be removed by Borel transformation .For these reasons, in order to obtain physical results the 
authors of ref. [^^^ considered various values of /i^ in the interval 0.1 < /x^ < 1 GeV^ and included 
uncertainties arising from this procedure into the error. From the presented above point of view 
such an approach is not legitimate. In this case p(si, s^) is proportional to Si, S2'- 

where 6 is a constant. In the model of hadronic spectrum accepted in Sec. 2, we have after the 
Borel transformation and using eq. (53) 
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r(M2) = 2P / dS2 / rfsi 

JO io 



Jo 



si - S2 



(87) 



Eq.(27) essentially differs from the corresponding expression for the bare loop contribution in 
ref.l^^': e.g., the integrand in (87) is positive, while in t^^l it is negative in the main region of 
integration. Of course, the QCD sum rule calculation in this case has a serious drawback: the 
continuum threshold W"^ dependence of the result is not in the form of a small correction of the 
type exp{—W'^/M'^) -C 1 at W"^ ^ M^, but much more strong. This is a direct consequence of 
high (equal to 5) dimension of the operators f/^, V^v^^- It is clear that the higher is the dimension 
of the considered operator, the stronger will be dependence on the continuum threshold and less 
certain the results of the QCD sum rules calculations. It must be emphasized that for operators 
of high dimensions the loop diagrams are in principle nonrenormalizable, the role of excited 
states in the physical spectrum increases and the determination of the lowest state contribution 
becomes impossible. 

Eq.(87) must be taken instead of the contribution of the bare loop diagram, used in '^'^1. A 
similar procedure must be also applied in the case of other terms in the sum rules I^^l containing 
an ultraviolet cut-off. 

BBK accounted for operators up to dimension 8. With the corrections described above, the 
sum rules found by BBK have the form 

((C/^^)) + A^' M' = e-'/^' I ^ - r 'ds e-^/^' + sin—) 





M'a'ln^ + ^TT^ nM^ - \ (88) 

((yivs^) + = e-'/M^ J _^ ^ Z"^' sds e"^/^' + sin—) 

"^^^ 135 n Jo s 

-- bM^E, [l—]-^^ M^a^ { In^ + O.d] + -vr^ RM^ - -my] (89) 
9 / 27 TT V M2 ; 27 9 ° J ^ ' 

in the notation(57),(58),(73). Here 11 and R are unduced by the external field v.e.v.'s which 
were estimated by BBK as 11 3.10"^ GeV^, R ~ 1.10~^ GeV^. Using these sum rules, one gets 

((f/^^)) ^ 0.14 GeV^ (90) 

((\/^^)) ^ 0.254 GeV^ (91) 

The substitution of (90), (91) into (80) gives the value of twist 4 correction to the Bjorken sum 
rule presented above in (10). The examination of the sum rules (88), (89) shows, however, that 
the final result (90) comes almost entirely from the contribution of the last term in OPE in (88) 
- the operator of dimension 8. Recently, A.Oganesian t^^l had calculated the contribution of 
factorizable v.e.v's of the operators of dimension 10, < | | >< qq \ 
g'^ < \ q G^^, X^q | >2 to the sum rule (88) and had found that by absolute value it is 
equal to (90) but has opposite sign. Of course, we cannot believe in this statement either: it 
means only that the results of the calculations are unstable and the value (10) characteizes the 
answer by the order of magnitude only. 
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4 Gerasimov, Drell-Hearn (GDH) Sum Rules. The Esti- 
mate of Higher Twist Corrections to F^^ Using Inter- 
polation between Vp^n and GDH Sum Rules. 

The real photon-nucleon forward scattering amplitude with nucleon spin flip is expressed through 
one structure function. In lab. system we can write 

e?\Tik)spinfHpek^ = i ^ ^iki ei^^ siS^{y,^) (92) 

777- 

where eS^^ and e*^^^ are polarizations of initial and final photons. At high energies, z/ — oo 
according to Regge theory the behaviour of 0) is determined by the exchange of the a\- 
Regge pole t^*^]; 

-Si(z/,0).^oo~^^""^^°^-' (93) 
Since Q;ai(0) ~ —0.3 — 0.0 the unsubtracted dispersion relation can be written for S\{y^ 0) 

S,{u) = J ^.p^^ (94) 



where Gi{i',Q'^) is the spin-dependent structure function. Consider the limit i/ ^ in (94). 
According to the F.Low theorem the terms proportional to and in the expansion in powers 
of u of the photon-nucleon scattering amplitude at small u are expressed via static characteristics 
of nucleon, its charge and anomalous magnetic moment 
The calculation gives 

Si{iy),-.o = (95) 

where k is the nucleon anomalous magnetic moment: Kp — 1.79, Kn — —1.91. 
Prom (94), (95) the GDH sum rule follows: 

oo 

I^G^iu,0)^-\K' (96) 



Till now no direct check of GDH was done. Only indirect check of (96)was performed, where in 
the l.h.s of (96) the parameters of resonances, obtained from the nN scattering phase analysis, 
where substituted. In this way with resonances up to 1.8 GeV it was obtained I^^'^o] (cf.also I^^l) 

l.h.s of (96) r.h.s. of (96) 

proton - 1.03 - 0.803 (97) 

neutron - 0.83 - 0.913 

The l.h.s. and the r.h.s. of (96) are not in a good agreement - a nonresonant contribution is 
needed. The direct check of the GDH sum rule would be very desirable! 

An important remark: the forward spin dependent photon-nucleon scattering amplitude has 
no nucleon pole. This means that there is no nucleon contribution in the l.h.s. of GDH sum rule 
- all contributions come from excited states: the GDH sum rule is very nontrivial. 

In order to connect the GDH sum rule with Tp^n{Q^) consider the integrals '^^^ 

/■oo fji/ 

Ip,n{Q')=L —Gl,p,n{^,Q') (98) 
J (D /2 
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It is easy to see that at large Q"^ 

IpAQ') = ^ ^pAQ') (99) 

and at = (98) reduces to the GDH sum rule. The dependence of Ip,n{Q^ is plotted 
in Fig. 11. 



Fig.ll. The connection o the GDH sum rules with the sum rules at high 
- qualitative dependence of /p,„ and Ip — In- 



The case of Ip is especially interesting: Ip is negative at = and positive at large Q^, 
what indicates to large nonperturbative corrections. In l^^l the VDM based interpolation model 
was suggested, describing Ip,n{Q^) in the whole domain of Q^. The model was improved in l^^l, 
where the contributions of baryonic resonances up to 1^ = l.SGeV, taken from experiment, 
where accounted. The model has no free parameters, besides the vector meson mass, for which 
the value iiy = 0-6GeV^ was chosen. Using this model it is possible to calculate the higher twist 
contributions in (2), (12). The results are presented in Table 1, as the ratio of asymptotic F"* 
with power corrections excluded to the experimentally measurable F at given Q'^ l^^l . 

Table 1. 

Higher twist corrections in GDH sum rule + VDM inspired model. 



Q^GeV^) 


2 


3 


5 


10 




1.44 


1.29 


1.18 


1.08 




1.30 


1.20 


1.13 


1.06 


-pas /p 


1.45 


1.29 


1.18 


1.08 


pos /p 

^ p+n/ ^ P+n 


1.47 


1.31 


1.19 


1.08 



The power corrections, given in Table 1 are essentially larger (except for the case of neutron), 
than the values (10), (86) found in l^^l. It must be mentioned that the accuracy of the model 
in the domain of intermediate Q^, where it is exploited, is not completely certain. So, 1 will 
consider in what follows the values of the power corrections presented in Table I as a limiting 
case of large higher twist corrections. 
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5 Comparison with Experiment. 



When comparing the sum rules (2), (12) with experiment I consider two hmiting variants of 
perturbativc corrections: small with A3 = 200MeV and large with A3 = 400MeV. (a^ is computed 
in 2- loop approximation, it is assumed that the number of flavours Nf = 3). For higher twist 
correction I also consider two limiting options: small (S), given by (10), (86) and large (L), 
determined by the data of Table 1. The contribution of gluons Ag{Q'^) in (12) will be found in 
the following way. Let us assume, that at 1 GcV the quark model is valid and ^^(IGeV^) = 
in eq.(18). Taking S=0.3, what is a reasonable average of the data, we have from (18) 



Ag{lGeV'^) = 0.35 
The dependence of Ag can be found from the evolution equation l^^l 



(100) 



1 + 



2N^ 
bn 



Q;,(Q^)-a,(/x2)]|A^(/x2)+ 



+ 



4 



- 1 



S(/.^), 



(101) 



where 5=11- (2/3)A^/ = 9, /i^ = IGeV^ and i:{lGeV^) ^ 0.3. As the calculation shows, the 
change of scale at which the quark model is assumed to work (say O.bGeV^ instead of IGeV^) 
or the of use shghtly different S(/i^) in (101) only weakly influence the results for E and As, 
obtained from experimental data. 

I consider the following experimental data (Table 2). 



Table 2. 

The experimental data on Fp^, 



Experimental 
group 


The target 


F(p, n) 


Mean Q^{GeV^) 


EMC [63] 


P 


Fp = 0.126 ±0.010 ±0.015 


10.7 


SMC [64] 


P 


Fp = 0.136 ±0.011 ±0.011 


10.5 


E143 [65] 


P 


Fp = 0.127 ±0.004 ±0.010 


3 


E142 [66] 


He" 


F„ = -0.022 ±0.011 


2 


SMC [68] 


d 


Fd = 0.034 ±0.009 ±0.06 


10 






Fp + F„ = 0.073 ± 0.022 


10 


E143 [67] 


d 


Td = 0.042 ± 0.003 ± 0.004 


3 






Fp ± F„ = 0.0908 ± 0.006 ± 0.008 


3 



The last column of Table 2 gives the average Q"^ in each experiment. In comparison with 
experiment the perturbativc and higher twist corrections, as well as Ag{Q'^) contributions are 
calculated for these Q'^. Experimentally, are different in different x-bins (higher at larger 
x). This effect is not accounted in the calculation. The ratio of term to term in perturbativc 
corrections is of order of 1 at A3 = 400Mey and Q'^ ^ 2 — SGeV'^ (as well as aj/a^ estimate). 
For this reason we introduce in these cases an additional error equal to the a^. The values of E 
and As calculated from comparison of experimental data with eq.l2 are shown in Table 3. (The 
errors are summed in quadrature). 
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Table 3. 

Determination of E and As from experimental data 



Experiment 


A3 


High 


E 


As 




MeV 


twist 








200 


S 


0.21 ±0.17 


-0.13 ±0.06 


EMC 


400 


S 


0.29 ±0.17 


-0.10 ±0.06 


P 


200 


L 


0.285 ±0.17 


-0.10 ±0.06 




400 


L 


0.37 ±0.17 


-0.07 ±0.06 




200 


S 


0.30 ±0.14 


-0.10 ±0.05 


SMC 


400 


S 


0.39 ±0.14 


-0.07 ±0.05 


P 


200 


L 


0.39 ±0.14 


-0.07 ±0.05 




400 


L 


0.47 ±0.14 


-0.04 ±0.05 




200 


S 


0.28 ±0.10 


-0.10 ±0.03 


E143 


400 


s 


0.42 ±0.10 


-0.06 ±0.03 


P 


200 


L 


0.57 ±0.10 


-0.006 ±0.03 




400 


L 


0.71 ±0.10 


0.04 ±0.03 


E142 


200 


s 


0.60 ±0.12 


0.003 ±0.04 


n 


400 


s 


0.57 ±0.12 


-0.005 ±0.04 




200 


L 


0.64 ±0.12 


0.016 ±0.04 




400 


L 


0.61 ±0.12 


0.008 ± 0.04 




200 


S 


0.27 ±0.10 


-0.11 ±0.03 


SMC 


400 


S 


0.33 ±0.10 


-0.09 ±0.03 


d 


200 


L 


0.29 ±0.10 


-0.10 ±0.03 




400 


L 


0.34 ±0.10 


-0.08 ±0.03 




200 


S 


0.37 ± 0.0(3 


-0.07 ± 0.02 


E143 


400 


s 


0.44 ±0.06 


-0.05 ±0.02 


d 


200 


L 


0.48 ±0.06 


-0.04 ±0.02 




400 


L 


0.54 ±0.06 


-0.015 ±0.02 



Remark: the contribution to E of the term proportional to Ag is approximately equal to 0.06 in 
the case of A3 = 200MeV and 0.11 in the case of A3 = AOOMeV. 

If we assume that all the analysed above experiments are correct in the limits of their quoted 
errors (or, may be, 1.5 st. deviations), then requiring for the results for E and As from various 
experiments to be consistent, we may reject some theoretical possibilities. A look at the Table 3 
shows that the variant A3 = 400MeV", L (a contradiction of E143, p and SMC, d results for E) 
and, less certain, the variant A3 = 200MeV, S (a contradiction of E142, n and SMC, d) may be 
rejected. 

Consider now the Bjorkcn sum rule. For comparison with theory I choose combinations of the 
SMC data - on proton and deutcron, the E143 data - on proton and dcutcron and the E143 data 
on proton and the E142 on neutron i^He). The results of the comparison of the experimental 
data with the theory are given in Table 4. 
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Table 4. 

Comparison of the experimental data witli tlie Bjorken sum rule 



Combination 


(Tp ^n) exper. 


A3{MeV) 


High 


— ^n)th 


of experiments 






twist 








200 


S 


0.193 


SMC, p 


0.199 ±0.038 


400 


S 


0.186 


SMC, d 




200 


L 


0.180 






400 


L 


0.173 


E143,p 




200 


S 


0.182 




0.163 ±0.010 ±0.016 


400 


s 


0.168 


E143,d 




200 


L 


0.145 






400 


L 


0.134 


E143,p 




200 


S 


0.182 


E142 n{^He) 


0.147 ±0.015 


400 


S 


0.168 


recalculated to 




200 


L 


0.145 


Q2 = 3GeV^ 




400 


L 


0.134 



From Table 4 we see again some indications for rejection of variants A3 = 400 MeV, L and, less 
certain, A3 = 200MeV, S.(In the first case there is a contradiction of the theory with the E143,p 
and d data, in the second - with the E143,p , E142, n data). 

At existing experimental accuracy it is impossible to choose from the data the true values of 
A3 and twist-4 correction . My personal preference is to the variant A3 = 200MeV and to the 
value of twist-4 corrections 3 times smaller than given by the GDH sum rule + VDM inspired 
model and, correspondingly, bp = 0.04 in (12), i.e., 2.2 times larger than the BBK result. The 
argument in the favour of such choice is that at larger A3 there will arise many contradictions with 
the description of hadronic properties in the framework of the QCD sum rules. The recent SLAC 
data [^^1 on gi{x, Q"^) Q^-dcpcndencc indicate that r^^/Fp = 1 + Cp/Q^, Cp = 0.25 ± 0.15 what is 
compatible with the estimate above. In this case all experimental data except for E142,n, are in 
a good agreement with one another and the values of S and As averaged over all experiments, 
except for E142,n are 

E = 0.35 ±0.05 As = -0.08 ± 0.02 (102) 

(see also l™' where the values close to (102) were obtained). 

The values of E and As obtained from the E142,n experiment at such a choice of A3 and 
twist-4 corrections, are different: 

E = 0.61 ±0.12 As = 0.01 ±0.04 (103) 

It is impossible to compete (102), (103) by any choice of A3 and higher twist correction. Perhaps, 
this difference is caused by inaccountcd systematic errors in the E142 experiment. 

A remarkable feature of the result (102) (as well as of the data in Table 3) is the large value 
of I As I - the part of the proton spin projection carried by strange quarks. This value may be 
compared with the part of the proton momentum carried by strange quarks 

^ I dxx [s+{x) + s_{x)] = 0.026 ± 0.006 t^^l, 0.040 ± 0.005 l^^^ (104) 
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The much larger value of | As | in comparison with V2 contradicts the standard parametriza- 
tion 

s+{x) + = A x-'^ (1 - xf 

s+(x) - s^x) ^ B x"^ (1 - x)^ (105) 

and requirement of positiveness of s+ and s_, if a ~ l(pomeron intercept) and 7 < (oi 
intercept). Large | As | and small Vf means that the transitions ss ^ uu + dd are allowed 
in the case of the operator j^s and are suppressed in the case of the quark energy-momentum 
tensor operator Q^^, corresponding to the matrix element V2. Such situation can be due to 
nonperturbative effects and the instanton mechanism for its explanation was sugg ested [^3]. It 
must be mentioned that in the suggested by Brodsky '^^1 more refined parametrization, which 
takes into account the fact that at x — > 1 g_(x)/g+(x) ~ (1 — x)^, the contradiction weakens. 
Improvement of experimental accuracy is necessary in order to be sure that the inequality 
I As 1^ V2 indeed takes place. 



6 The calculations of the polarized structure functions 
by QCD sum rule method 

1 recall the basing points of the calculation of nucleon structure functions in the QCD sum rule 
approach (for details see '^^1. Consider four-point correlator 

T^^{p,q) ^-i J d^xd^yd^z e'«^ e^^'^^-^^x 

X < I T {77(y), jj(x), i±(0), f]{z)] I > (106) 

where 

and r) is given by (36). (In order to separate u- and d quark distributions it is convenient to 
consider proton scattering). Suppose that < 0,p^ < 0,g^ = —Q^, is large enough, 
~ 10 GeV'^, >>| I and retain only the leading terms in the expansion over /Q"^. (This 
corresponds to account of only twist 2 contributions). Assume also that | |>> R~^i where Rc 
is the confinement radius and perform OPE in for the discontinuity of T^^, in the s-channel 



Tt,ij>\ 9', s + k)- T±(p', q\ s - it) (107) 



At the first sight it seems, that OPE is not legitimate here because T^j, is the forward scattering 
amplitude in which large distances in the t-channel are of importance. This, indeed, would be 
the case if we would consider the exclusive object Re T^i,. But for Im T^i, which in fact is an 
inclusive observable, the situation is completely different. Consider first the bare loop diagram 
for Im T^^ (the contribution of dimension operator in OPE) - Fig. 12. 
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Fig. 12 The bare loop diagram for Im T^^. The crosses mean on mass shell propagators. 

The direct calculation of the diagram gives that the virtuality of the active quark on which the 
scattering proceeds is equal to {x — Q^/2z/, v = pq) 

k^^p^x-^ (108) 
1 — X 

and k']_ is of order fc^ ~| | a;(l — x). Therefore, k'^ ~ p'^x. (Strictly speaking, this statement 
refers only to terms in the amplitude, singular in p^, but just these terms are of interest for 
us since the Borel transformation in killing the regular in terms will be performed later). 
Since it was assumed that | p^ |>> we come to a conclusion that at not small x,\ k"^ \ and 
k^ are large and OPE is legitimate. 

A general proof of the above mentioned statement follows from the fact that at large | | , 
I I the nearest to zero singularity in t of the function Im T{p^, q^, s, t) is determined by the 
boundary of the Mandelstam spectral function which is found to be 

S 1 — X 

So, at t = 0, large | p^ \ and not small x, we are far from the boundary of the spectral function 
in t-channel and OPE is valid. 

The statement that the method does not work at small x is evident beforehand: this is the 
Regge domain, where OPE cannot work. 

The estimate of active quark virtualities (108) have much more general meaning, beyond QCD 
sum rule appoach. In lepton-hadron scattering | p^ | is of order of QCD scale, | p^ |~ 1 GeV^. 
Then at small x \k'^ \ are small and we are in the nonperturbative domain of QCD. This means 
that the interaction of quarks with nonperturbative vacuum fields, especially gluonic fields, are 
important. The direct calculation in QCD sum rule approach supports this expectation. 

The method in view is also invalid at x close tol, 1 — a;<<l. This is evident, because it is 
a resonance region. Finally, we restrict ourselves to intermediate ~ 5 — 10 GeV^, since the 
evolution of the structure functions will not be accounted. 

The calculation proceeds in the standard way of the QCD sum rule method. In the QCD 
side of the sum rule Im is calculated by OPE with the account of v.e.v.of various operators. 
For the nonpolarized case the gluonic condensate and the term ~ < | | >^ were 
accounted For polarized structure functions gi and (72, besides the bare loop, only the 

contribution of the term ~ < | gg | >^ was calculated since it is expected that at 
intermediate x, where the results are correct, the contibution of gluonic condensate is small. 
The v.e.v.< | gg | >^ do not contribute to the sum rule, because it is concentrated x — 1 
(proportional to 5(1 — x) in twist 2 terms) - outside of the applicability domain of the method. 
(Only chirality conserving structure were considered in I'''^''''^]). The hadronic side of the sum rule 
is represented schematically by Fig. 13. 
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Fig. 13. The representation of hadronic part of the sum for determination of the nucleon structure functions: a) 

the double nucleon pole term which gives the desired structure function; b) the background term, when in the 
course of the deep inelastic scattering N undergoes transition to excited state N*; c) continuum. 



The Borcl transformation is appUed to both sides of the sum rule. The continuum contribution 
is suppressed by the Borel transformation, it is approximated by the bare loop (Fig. 12) and 
transferred to the QCD side. The background term - Fig. 13b - in its dependence on the Borel 
parameter differs from the nucleon term by an additional factor M^.This circumstance 
permits one to kill this term in the same way, as it was killed the constant A in (56). 

In the case of polarized structure functions gi and g2 it was shown t'^'^l that for the bare 
quark loop the Bjorken and Burkhardt-Cottingham '^^1 sum rules are fulilled. It was found 
that for the function gi{x) the results are reliable in a rather narrow domain of intermediate 
x: 0.5 < X < 0.7. In this domain the contribution of li-quarks g^ is much larger than d-quarks, 
91 ^ gf- Therefore, gi (4/9) g'^. g'^ was calculated and for the mean value of gi in this 
interval it was obtained (at ~ 5 — lOGeF^): 

^i(0.5 < X < 0.7) =0.05±50% (110) 

The large uncertainty in (110) results from the large contribution of nonleading term in OPE 
and from large background at the phenomenological side of the QCD sum rule. The E143 proton 
[^^1 and deuteron [^^1 data (the latter under assumption that g'^ is small in this interval of x) give 
roughly the same values: 

^1 (0.5 < X < 0.7) = 0.08 ± 0.02 (111) 
This value is in a good agreement with the SMC result l^^l 

^i(0.4 <x< 0.7) = 0.08 ± 0.02 ± 0.01 (112) 

and with theoretical expectation (110). 

For the case of the structure function g2 only gl^ can be calculated at 0.5 < a; < 0.8, the 
calculation of g2 fails because of large contribution of nonleading terms in OPE. If we assume 
that hke in the case oi gi, | 5'2 l"^! I; ^^en ['''^1 

^2(0.5 <x< 0.8) = -0.05 ± 50% (113) 

The E143 data ''''^1 in this interval of x are: 

52(0.5 <x< 0.8) = -0.037 ± 0.020 ± 0.003 (114) 

in a good agreement with (113) 

The serious disadvantages of the QCD sum rules calculations of the structure functions are: 
1) large contribution of nonleading terms in OPE; 
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2) large background terms from inelastic transitions N — > N* (Fig. 13b) at the physical side 
of the sum rules. In order to kill these terms we are forced to differentiate the sum rules over 
the Borel parameter M^. This operation increases the role of nonleading terms in OPE and 
continuum and deteriorates the accuracy of the sum rule. A possible way to overcome this 
drawback of the method is to start from nonforward scattering amplitudes qi ^ q2,Pi 7^ P2 
to use the double Borel transformation in pi and p"^. The calculations in such approach will be 
much more complicated, but, may be the game is worth of candles. 

7 Calculation of Chirality Violating Structure Function 

hi{x) by QCD Sum Rules 

As is well known, all structure functions of the twist two - Fi{x), F2{x), gi{x), which are 
measured in the deep-inelastic lepton-nucleon scattering, conserve chirality. Ralston and Soper 
first demonstrated that besides these structure functions, there exists the twist-two chirality 
violating nucleon structure function hi{x). This structure function does not manifest itself in 
the deep inelastic lepton-hadron scattering, but can be measured in the Drell-Yan process with 
both beam and target transversally polarized. The reason of this circumstance is the following. 
The cross section of the deep inelastic electron(muon)-hadron scattering is proportional to the 
imaginary part of the forward virtual photon-hadron scattering amplitude. 



Fig.l4 

a) Deep inelastic lepton-hadron scattering, the 
quark chiralities are conserved. Solid lines are 
quarks, wavy lines are virtual photons, R{L) de- 
note right (left) chirality of quarks; b) Drell-Yan 
process with chirality of quarks flipped. 

At high photon virtuality the quark Compton amplitude dominates, where the photon is 
absorbed and emitted by the same quark (Fig. 14a) and the conservation of chirality is evident. 
The cross section of the Drell-Yan process can be represented as an imaginary part of the 
diagram. Fig. 14b. Here virtual photons interact with different quarks and it is possible, as is 
shown in Fig. 14b, that chirality violating amplitude in Drell-Yan processes is not suppressed at 
high Q'^ in comparison with chirality conserving ones, and consequently, corresponds to twist 
two. This amplitude, corresponding to target spin flip, has no parton interpretation in terms 
of quark distritutions in the helicity basis and, as was shown by Jaffe and Ji 1^°^ can be only 
represented as an element of the quark-quark density matrix in this basis. 

However, hi{x) can be interpeted as a difference of quark densities with the eigenvalues 
-1-1/2 and —1/2 of the transverse Paili-Lubanski spin operator sj^75 in the transversely polarized 
proton. It this basis hi{x) can be described in terms of standard parton language. 

The proton structure function hi{x) can be defined in the light cone formalism as follows 

i I f^e^^"" <P,s\ '0(O)(T^^75V'(An) | p, s >= 2[hi{x, Q^){s^f,p^ - s±^p^)+ 
+hL(x, Q'^)rin?{pf,n^ - p^nf,){sn) /i3(x, (5^)m^(s±^n^ - sx_^n^)\ . 

Here n is a light cone vector of dimension (mass)~^, = 0,n+ = 0,pn = l,p and s are the proton 
momentum and spin vectors, p^ — m^, — —l,ps — and s — {sn)p+ {sp)n + s±, hL{x, Q^) and 
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h3{x, Q"^) are twist-3 and 4 structure functions. For comparison in the same light cone notation 
the standard structure function Q^) is given by 



r d\ 

I 2^"^'^" ^ ^' ^ ' VS(0)7MV'(An) I p, s >= 4[Fi(x, Q'')p^ + MV4(x, Q')nJ (116) 

(Eqs. (115), (116) are written for one flavour). 

Basing on the definitions (115), (116) an inequality was proved in ^° 

q{x) > hl{x) , (117) 

which holds for each fiavour q = u,d,s. (Here hi{x) is the fiavour q contribution to hi{x).) 
Recently, Soffer ^^^^ (see also t^^l derived an inequality 

\K\<[u{x)+g'^{x)]/2 (118) 

hi{x) can be also represented through a T-product of currents '^^1 

T^{p,q,s) = i J d^xe^'i^ <p,s\ {l/2)T{j^,{x), j{0) + j{x), Ji,,{0)} \p,s> , (119) 

where j^i5{x) and j{x) are axial and scalar currents. 
The general form of Tn{p, q, s) is 

T,{p, q, s) = (^s, - ^g, j h{x, Q') + [p, - ^ j {qs)h{x, Q') + 

+ £tiv\aPvq\Sa{qs)l2{x,q^) (120) 
(only spin-dependent terms are retained). It can be proved that 

hi{x,Q'^) = - -7m/ii(x,Q^) . (121) 

TT 

As is clear from (115) or (119) hi{x) indeed violates chirality. 

Since hi{x) violates chirality one may expect that in QCD it can be expressed in terms of 
chirality violating fundamental parameters of the theory, the simplest of which (of the lowest 
dimension) is the quark condensate. Basing on this idea hi{x) calculation was performed '^^1. 
The main difference in comparison with the structure functions considered in Sec. 6 is that in the 
case of hi{x) determination the chirality violating structure is studied. As a result, hi{x) was 
found to be proportional to the quark condensate < | mm | > with the correction term in OPE 
proportional to the mixed quark-gluon condensate < | ua^^, G^^ X^u | >. It was obtained 
that for proton h^{x) » hf{x) and as a consequence hi{x) ^ {4:/9)hi{x). The calculations in l^^l 
are valid at 0.3 < x < 0.6. The extrapolation in the region of small x can be performed using the 
Regge behaviour hi{x) ~ the extrapolation in the region o large x, using the inequalities 

(117), (118). The final result of the calculation of hi{x) with the extrapolation is shown in Fig. 15. 
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Fig. 15 The w-quark contribution to the proton structure 
function hi (x) based on the QCD sum rule calculation 
[^^1 at intermediate x, 0.3 < x < 0.6. At x < 0.3 
an extrapolation according to Regge behaviour was performed. 
The dashed line represents the Soffer inequality (118). 

It is expected that the accuracy of hi{x) determination is about 30% at a; = 0.4 and about 50% 
Bit X — 0.6. The inequahty hi{x) > gi{x) suggested in was confirmed. Numerically, hi{x) is 
rather large, that gives a good chance for its experimental study. 



8 Conclusions 

The experimental study of the nucleon spin structure, where very impressive results were ob- 
tained, triggered many theoretical investigations. As result, we understand now much more 
about the internal content of nucleon and, even generally, about the structure of QCD. We 
know, that EUis-Jaffe sum rule is not the last word in the problem of the nucleon spin con- 
tent: gluons and strange quarks are of importance in this problem. The connection of gluon 
contribution to the nucleon spin to the anomaly is clarified. More clear becomes the role of 
nonperturbative phenomena in QCD. In this aspect there is a very important indication that 
the part of the proton spin projection carried by strange quarks | As | is much larger than the 
part of the proton momentum V2 carried by strange quarks, | As |3> Vf^. (A confirmation of 
this with a better accuracy is necessary!). If confirmed, this statement indicates a nontrivial 
dynamics of QCD vacuum and its explanation is a challenge for a theory. It would be very 
desirable if experimentally it would be possible in the near future: 

1. To increase the accuracy by 2-3 times. 

2. To study the Q^-dependence of gi{x, (5^)(in separate bins in x). 

3. To go to higher Q"^: probably this can be done at HERA. (The experiment at higher 
will be informative if only its accuracy will not be worse than the existing ones). 

4. To have better data in the domain of small x - probably, this also can be done at HERA. 
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5. To perform measurements of two-jets events in polarized deep-inelastic scattering, which 
will give information about /S.g. 

6. To have more data on s-quark distribution in nonpolarized nucleon, expecially at x < 0.1. 

7. To perform new experiments on elastic z/p-scattering from which one can find the combi- 
nation Am — A(i — As, as well as elastic ep scattering with separating the Z-exchange term. 

8. To have better data on g2. 

9. To measure hi{x). 

10. To perform a direct check of the GDH sum rule. 

The problem confronting the theory, which, 1 believe, could be solved in the near future, are: 

1. The lattice calculation determination of the induced by external field v.e.v's, especially 
those which cannot be calculated by the QCD sum rule approach. 

2. The study of ss and uu + dd mixing in the nucleon (on lattice and in the instanton liquid 
model), especially for the case of the axial and energy- momentum tensor operators. 

3. Improvement of the calculation of the structure functions by the QCD sum rule method. 

4. To achieve a better understanding of nonperturbative effects in the structure functions at 
small (but not very small) x ~ 10~^. 
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